CONFORMAL INVARIANCE OF LOOPS IN THE 
DOUBLE-DIMER MODEL 

RICHARD KENYON 

Abstract. The dimer model is the study of random dimer covers 
(perfect matchings) of a graph. A double-dimer configuration on 
a graph Q is a union of two dimer covers of Q. We introduce 
quaternion weights in the dimer model and show how they can 
be used to study the homotopy classes (relative to a fixed set of 
faces) of loops in the double dimer model on a planar graph. As 
an application we prove that, in the scaling limit of the "uniform" 
double-dimer model on Z 2 (or on any other bipartite planar graph 
conformally approximating C), the loops are conformally invariant. 

As other applications we compute the exact distribution of the 
number of topologically nontrivial loops in the double-dimer model 
on a cylinder and the expected number of loops surrounding two 
faces of a planar graph. 



1. Introduction 

A dimer cover, or perfect matching, of a graph Q is a set of edges 
with the property that each vertex is the endpoint of exactly one edge. 
A double-dimer configuration is a union of two dimer covers, or 
equivalently a set of even-length simple loops and doubled edges with 
the property that every vertex is the endpoint of exactly two edges 
(which may be doubled). 

The uniform dimer model (the uniform measure on dimer covers) 
on the square grid Z 2 has been the subject of much research, starting 
with the exact enumeration of j5j [18] and culminating in the limit shape 
theorems in, successively, [TJ [TTJ [TU] , and conformal invariance results 
in [8]. It was conjectured in [T3] (see also [IE] ) that the loops in the 
double-dimer model also have conformally-invariant scaling limits, and 
are related to SLE4. We prove here the conformal invariance of the 
loops. 

If a graph Q has positive edge weights v: E — > IR>o, there is a natu- 
ral probability measure /i„ on dimer covers which gives a dimer cover 
a probability proportional to the product of its edge weights. From 
Hv ® we similarly get probability measures on double-dimer con- 
figurations. More general double-dimer probability measures can be 
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constructed by taking //^ <g> /i^ 2 for different z/i,z/ 2 , or even z/ 2 = \>\ 
when 1/1 is complex- valued (in this case \i Vx <8> //^ may be a probability 
measure on double-dimer configurations, thought of as collections of 
loops, even though \i Vx and not probability measures). 

The main innovation in this paper is the introduction of quater- 
nionic (or more generally SL 2 (C)) weights for the double-dimer model, 
a tool which is not available for the single-dimer model. The Kasteleyn 
theorem on evaluating the weighted sum of dimer covers via determi- 
nants can be extended to apply in this case; one must use the so-called 
quaternion-determinant, or g-determinant [31 [H] of a self-dual quater- 
nionic matrix. The g-determinant of the quaternionic Kasteleyn matrix 
is then the partition function of the double-dimer model (Theorem [2] 
below). 

Another important tool in the proof of conformal invariance is a 
theorem of Fock and Goncharov [I], Theorem [3] below, connecting the 
set of simple closed curves (or more generally laminations) on surfaces 
with the SL 2 (C) representation theory of its fundamental group. 

On a multiply-connected planar domain U we show that, for a graph 
conformally approximating U (in the sense that discrete harmonic func- 
tions converge to harmonic functions on U), the distribution of the ho- 
motopy types of the homotopically nontrivial double dimer loops only 
depends on the conformal type of U. See Corollary [5] for the case when 
the graph is essentially U fleZ 2 , and Section [8] for the extension to other 
graphs. 

As a limiting case of the above we can take U to be a bounded 
simply connected region from which a finite number of points Z\, . . . , z m 
have been removed; in this case we show that the distribution of the 
homotopy types of the nonperipheral, homotopically nontrivial loops 
in the double dimer model only depends on the conformal type of U. 
(A loop is peripheral if it is isotopic to a small loop surrounding exactly 
one of the punctures Zj.) See Section [7J 

We give two other applications of these ideas: a computation of 
the exact distribution of the number of topologically non-trivial loops 
in a double-dimer configuration on the square-grid on a cylinder (see 
Section|9|), and the expected number of double-dimer loops surrounding 



simultaneously two points z±,Z2 (see Section 10). 

In our main proof we approximate U with a discrete graph Q t C eZ 2 
with rather special boundary conditions called Temperleyan boundary 
conditions [7], (See Figure [2]). This issue of choice of boundary condi- 
tions is a delicate one. In fact, for naive choices of boundary conditions 
the conformal invariance will not hold, see [5]. There are nonetheless 
many other "good" choices of boundary conditions and it is possible to 
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show with the techniques in this paper that the conformal invariance 
holds for these graphs as well. Our choice in the current paper was de- 
signed to make the proof as simple as possible. The technical difficulty 
in the general case is showing convergence of K" 1 , the inverse Kaste- 
leyn operator. In the current case we can write K~ x in terms of the 
discrete Dirichlet Green's function whose convergence properties are 
well known. For other boundary conditions the corresponding discrete 
Green's functions are (usually) less well understood. 

2. Vector bundles and connections 

Let Q = (V, E) be a graph and W a vector space. A IV-bundle on 
Q is a vector space W v isomorphic to W associated to each vertex v of 
Q. The total space of the bundle is Wg = (B V £vW v ; a section is an 
element of Wg. A connection $ = {<f) e }e£E on a vector bundle is the 
choice of an isomorphism <p vv i : W v — > W v > for every edge e = vv', with 
the property that 0_ e = where — e is the edge e with the reverse 
orientation. The map <p vv i is the parallel transport of vectors at v to 
vectors at v'. 

Two connections $ and $' are gauge equivalent if there are iso- 
morphisms if) v : W v — > W v such that 4 , v 2 ( fi'v 1 v2 = < Aui«a'0vi for all adjacent 

Vi,V 2 . 

If 7 is an oriented closed path on Q starting at v, the monodromy 
of the connection along 7 is the product of the parallel transports 
along 7; it is an element of End(W v ). Gauge equivalent connections 
give conjugate monodromies; changing the starting point along 7 also 
conjugates the monodromy. 

If Q is embedded on a surface (in such a way that faces are con- 
tractible), a flat connection is a connection with trivial monodromy 
around faces, and thus around any contractible curve. If $ is a flat 
connection, given any contractible union of faces S, it is not hard to 
see that one can choose a connection gauge equivalent to $ which is 
locally trivial on S, that is, is the identity on each edge in S. 

It is convenient to represent a flat connection as in Figure |2] take 
a set of simple closed paths in the dual graph which support cocycles 
generating the cohomology of the surface (in other words, cutting the 
edges crossing these paths results in a contractible surface) and take 
flat connection supported on the edges crossing these paths. Any flat 
connection is equivalent to one of this form. We call such a dual path 
a zipper. 

In this paper we use only 1- or 2-dimensional complex vector bun- 
dles, with parallel transports in C* or SL 2 (C). In fact for the purposes 
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of studying the double-dimer measure it suffices to take unitary con- 
nections, with structure group Ui or SU2. 

3. DlMER MODEL 

3.1. The single-dimer model. 

3.1.1. Definition. Let Q be a graph and v : E — >■ IR>o a positive real 
weight function on the edges. Let M. = M.{Q) be the set of dimer cov- 
ers (perfect matchings) of Q. For m E M. define u(m) = rieem z/ ( e ) ^° 
be its weight. We define a probability measure \x on A4 where the prob- 
ability of m is proportional to v(m). The constant of proportionality 
is 1/Z where Z = Yl m eM z/ ( m ) i s called the partition function. 

Note that if v' is a different weight function, obtained from v by 
multiplying the edge weights at a given vertex by a nonzero constant A, 
then n is unchanged, since every dimer cover has weight multiplied by A. 
Compositions of such operations are called gauge transformations. 
Two weight functions are gauge equivalent if they are related by a 
gauge transformation. 

3.1.2. Line bundle interpretation. Now suppose Q is bipartite, with 
black vertices B and white vertices W. Define $ = & u , a connection 
on a line bundle on Q, by taking W = C 1 and e = v{e) if e is directed 
from black to white, that is, (f) e is multiplication by v[e) (we must then 
have 0_ e = z/(e) _1 ). If V\ and z/ 2 are gauge equivalent as weight func- 
tions, then their associated connections § Ul ,§ U2 are gauge equivalent 
as bundles. Thus fi depends only on the gauge equivalence class of 
gauge transformations do not change /z. 

3.1.3. Kasteleyn matrix. If Q is planar and bipartite, Kasteleyn showed 
[3] that the partition function Z is the determinant of the Kasteleyn 
matrix K = (K(w,b)) w£ y\; t b e B whose rows index white vertices and 
columns index black vertices, and K{w,b) = ±z/ w t = ±4> w b- Here the 
sign is chosen according to the "Kasteleyn sign condition" for bipartite 
graphs, that is, so that the number of minus signs around a face of 
length £ is f + 1 mod 2. We have Z = | det K\. 

In the case Q is planar but not bipartite, there is a generalization of 
this result also due to Kasteleyn: Z is the Pfaffian of an antisymmetric 
matrix K = (K(v , v')) v>v > e v (also called Kasteleyn matrix), indexed by 
all the vertices of Q, with K(v,v') = ±u(vv'). Here the sign is chosen 
by orienting the edges of Q in such a way that each face has an odd 
number of clockwise-oriented edges. 

3.2. Double dimer model. 
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3.2.1. Definition. Let Q be a general (not necessarily bipartite) graph 
and v\ , v 2 : E — > C two weight functions, not necessarily real. Given a 
pair (mi, 777.2) £ -M 2 , we associate a weight v(m\,m 2 ) = fi (777,1) ^2(^2) 
where Vi(m) are defined as above. The partition function is defined 

as 

(1) Z dd = Z{vx)Z{v 2 ) = v{m x ,m 2 ). 

(mi,m2)&M 2 

It is natural to group configurations according to the set of loops 
they form. Let Q = fl(G) be the set coverings of Q by collections of 
edges which form cycles of even length or doubled edges (so that each 
vertex is the endpoint of exactly two edges, which may be the same 
edge). We call Q the set of double dimer configurations on Q. To 
uGfiwe associate a weight 

v ( u ) = II vx{e)v 2 {e) Y[ i w i + w 2 ) 

doubled edges cycles 

where for a cycle 7 = (v% — >■ v 2 —)■••• — >■ — >■ ^i) we associate two 
weights 

^l(7) = ^l(^1^2)^2(w2^3)^l(^4)^2(^5) • • • 

and 

^2(7) = ^(viv^vxfav^v^vsv^v^v^) .... 
Lemma 1. Z dd = Z{v x )Z{v 2 ) = XLenK^)- 

Proof. A single element u € fl with k cycles corresponds to 2 k pairs 
(777,1, 7772) € A^ 2 : each cycle in u can be partitioned in two ways to get 
two dimer covers of it. □ 

If we are interested in constructing a probability measure on Q, we 
must (somehow) choose Ui, v 2 so that the weights of configurations are 
real and nonnegative. The reality can be guaranteed if we take for 
example v 2 — V\, since then w 2 = W\ for any cycle. Positivity can then 
be achieved if the arguments of the weights are sufficiently small. 

3.2.2. Computation. If Q is planar, from (JI]) we have 

(2) Z dd = PLfYiPf# 2 

where K\, K 2 are antisymmetric Kasteleyn matrices associated to v\, v 2 
respectively. 

If Q is bipartite and planar, we can instead use a single matrix K = 
(K(v,v')) VtV i e v indexed by all vertices, where K(w,b) = ±Vx(b,w) and 
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K(b, w) = ±1^2(10, b). The signs are again determined by the Kasteleyn 
condition for bipartite graphs mentioned above. In this case we have 

(3) Z dd = detK. 

Equations Q and ^ are special cases (for antidiagonal matrices and 
diagonal matrices, respectively) of the construction below. 

3.3. SL2(C)-bundle case. Let Q be a bipartite graph and v: E — > M>o 
a positive real weight function on the edges. Let <3> be an SL 2 (C)- 
connection on a C 2 -bundle on Q, that is, a connection where the parallel 
transports are in SL 2 (C). 

The weight of u G £1 is now defined to be 

(4) !/(«) = J] "(e) II Tr H, 

e cycles 

where the first product is over edges of u (and doubled edges are 
counted twice) and the second product is over nontrivial cycles (that 
is, not the doubled edges), w being the monodromy of the cycle in one 
direction or the other (note that Tr(w) = Tr(if -1 )). 

In order to associate a probability measure it is necessary that the 
trace of the monodromy on any loop which occurs in some u G Q be 
nonnegative. This can be arranged for example if the parallel trans- 
ports <p vv i are all sufficiently close to the identity. 

Associated to an SL 2 (C) connection on a bipartite planar graph is 
a Kasteleyn matrix K with entries in GL 2 (C), defined as follows. We 
define K(v,v') = K u (v,v')(f) vv i where K u is the usual weighted signed 
Kasteleyn matrix, and <p vv i are the parallel transports written in the 
standard basis for C 2 . 

Note that the matrix K is self-dual, that is K(v,v') = K(v',v)* 

where by * we mean the "g-conjugate": if A — ( a ^) then A* = 



d -b 
—c a 



Theorem 2. We have Z dd = XLen u i UJ ) = Qdetif. 

Here Qdet is the quaternion determinant of a self-dual matrix, de- 
fined by: 

(5) QdetK = ^ sgn(cr) J] ^Ti(K c ), 

trdSn cycles C 
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where the product is over cycles C of a, and Kc is the product of the 
elements of K along the cycle C, that is, if C is the cycle 

v i -> v 2 u„ -> Ui 

then .ft'c' — K VlV2 K V2 v 3 ■ ■ ■ K- VnVl . 

Proof. Each summand in (|5| for Qdet(i^) is zero unless a maps nearest 
neighbors to nearest neighbors, and thus is a double-dimer configura- 
tion; because K is bipartite all cycles must have even length. The 
weight of a nonzero term is equal to the product of its edge weights 
(with doubled edges counting twice) times the product of 1/2 the trace 
of the monodromy on all nontrivial (length bigger than 2) cycles. When 
we account for both orientations of each nontrivial cycle, its contribu- 
tion is the trace of its monodromy. 

It remains to verify that each nonzero term has the same sign. This 
is true by Kasteleyn's theorem (equation ^ above) for the trivial bun- 
dle, and thus (by for instance deforming continuously to the case of a 
nontrivial bundle) remains true in general. □ 

If K is an n x n self-dual matrix with entries in GL^C), let K be the 
2n x 2n matrix with entries in C, obtained by replacing each entry by 
its 2 x 2 block. Mehta [13] , see also pJj , proved that Qdet-fT is equal to 
the Pfaffian of the (antisymmetric) matrix ZK, where Z is the 2n x 2n 

matrix with blocks ( ^ n ) on the diagonal and zeros elsewhere. 



-1 / 

Since Q is bipartite, listing white vertices first, K has the form 
M N 



K = q J where M is the matrix with GL 2 (C)-entries with 

rows indexing white vertices and columns indexing black vertices. Thus 
Qdeti^ = det M where M is the n x n matrix obtained by replacing 
each entry in M by its 2 x 2 block. 

Later we will use the fact that the inverse of a self-dual matrix K 
of nonzero Qdet is well-defined, is both a left- and right-inverse, and is 
self-dual, see [3]. 

Note that if we choose 6 P — ( ^ ] where A = \/^44 and vie) 



" V ^( e ) 

vi(e)is 2 (e) then (J3J) is a special case of Theorem [2j 
The construction in this section also works for nonbipartite graphs 
if we can arrange that the trace of the monodromy around every odd- 
length loop is (then configurations with odd-length loops have zero 

weight). For example if we choose 4> e = ( ^ J where A = 
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and v(e) = \J v\[e)V2[e) then (2) is a special case of Theorem [2] for 
nonbipartite graphs. 

Question: Are there other, more interesting, connections in which the 
trace of the odd-length-loop monodromies are zero? 



4. Topology of paths 

An noncontractible simple closed curve on a surface E is a 
simple closed curve which does not bound a disk. A finite lamination 

on a surface S, possibly with boundary, is a finite set L = {71, . . . , 7^} 
of pairwise disjoint noncontractible simple closed curves 7i, considered 
up to isotopy. By this we mean that we consider two finite laminations 
to be equivalent if they are isotopic, that is, they have the same number 
of loops in the same homotopy classes. 

Let (Q, v, E) be a triple consisting of a bipartite graph Q, with edge 
weights u, embedded on a surface E. The loops in a double-dimer cover 
*jj G may or may not be contractible. We associate to u the finite 
lamination A (a;) which is the (possibly empty) set of noncontractible 
cycles, that is, those cycles which do not bound disks in E. 

Let no be the probability measure on double-dimer configurations fl 
defined for uj G fl by 

eSu; 

where Z is a normalizing constant, k is the number of nontrivial cycles 
(contractible or not) and again doubled edges are counted twice in the 
product (but do not count as cycles). This is the weight Q associated 
with the trivial connection. We refer to //q as the natural measure on 
double dimer configurations (since it arise from the uniform measure 
on single-dimer configurations). We wish to study the /io-probability 
of any finite lamination on E, that is, the probability that a /xo-random 
configuration u has X(u) of a specified type. 

Recall that a flat connection on (Q, E) is a connection whose mon- 
odromy around every contractible loop is the identity. The set of flat 
SL 2 (C)-connections modulo gauge transformations has the structure 
of an algebraic variety X = AT(E), isomorphic to the representation 
variety of homomorphisms of ^(E) into SL 2 (C), modulo conjugation. 

Let $ = {4>e} be a flat connection. By Theorem |2j the partition 
function Zdd = Zdd{&), which is a polynomial function of the entries of 
the e , is a nonnegative linear combination of terms ri 7 eL Tr(u> 7 ) where 
L runs through all possible finite laminations of (Q, E). In particular 
it is a function on X. 
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Let V = V(Q, S) be the vector space of functions of $ = {<p e } which 
are polynomial in the matrix entries and invariant under gauge trans- 
formations. This is a vector space of infinite dimension but finite in 
each degree. Fock and Goncharov proved 

Theorem 3 (|4J). Suppose £ has nonempty boundary. The elements 
Yl-y^L ^t)) as e/ements of V, are linearly independent, and in fact 
form a basis for V as L runs over all finite laminations. Moreover 
there is a natural inner product on V, coming from Haar measure on 
SL 2 (C), under which the functions Y\ 1&L Tr(w 7 ) are orthogonal. 

In particular one can extract, via an integral over the representa- 
tion variety X(Ti), the coefficient in Z^d of any given finite lamination. 
When the surface has nonempty boundary, 7i"i(E) is a free group and so 
X{T) is a just a product of copies of SL 2 (C). The corresponding inte- 
grals are integrals against spherical harmonics times the Haar measure 
on (SU 2 ) n . 

5. Scaling limit of the double-dimer path 

5.1. Scaling limit. Let U C C be a bounded multiply connected do- 
main with boundary consisting of a finite number of piecewise smooth 
curves Co, . . . , C m , with Co being the outer curve. Let Zi G Cj be fixed 
points, one on each boundary component. We assume (to simplify the 
proof of convergence of the Green's function) that the boundary of C 
is flat and horizontal in a neighborhood of each z^. Let 7$ be a simple 
path in the interior of U from z to Z{. We suppose these ji are pairwise 
disjoint except at z . 

For e > let U e be the graph with vertices U PI eZ 2 and edges 
connecting points at distance e. For each j > let ej be an edge of 
U e near Zj (that is, converging to Zj as e — > 0), whose vertices lie on 
the face containing Cj. Let xq be a vertex near Zq. See Figure [T] for an 
example. 

Let Q e be a bipartite graph with a black vertex for each vertex (ex- 
cept xq) and square face of U e (faces do not include the complementary 
components), and a white vertex for every edge except for the edges ej. 
Edges in Q e connect nearest neighbors, that is, "edge" vertices to "ver- 
tex" vertices and "edge" vertices to "face" vertices if the corresponding 
elements are adjacent. See Figure [2j The graph Q e is bipartite and has 
dimer covers (these are in fact in bijection with spanning trees of U e 
rooted at x with a certain property: their dual trees rooted on the 
outer face have branches from the dual vertices in the Cj which lead in 
the direction of the e,; see [T2|. 
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Figure 1. The region U and graph U e with marked 
edges e, and vertex xq. 
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Figure 2. The associated graph G e . The "zippers" are 
the thicker edges (in blue). 



Let fi e be the natural measure on double-dimer configurations Q(G e ) 
of Q e , that is, in which a configuration with c loops has a probability 
proportional to 2 C . 

For a double-dimer configuration u G Vt(Q t ) let L = L(u) be the 
finite lamination of U consisting of the noncontractible loops of u. 

We prove that the /^-distribution of L converges as e — > and 
the limit only depends on the conformal type of the marked surface 
(U, {z , . . .,z m }). 

5.2. Conformal invariance. Let $ be a flat SL2(C)-connection on 
{Qe, U). Let Ai G SL 2 (C) be the monodromy of a path in Q e isotopic to 
one traversing 7^, then running counterclockwise around the boundary 
of Ci, and then traversing 7^ back towards z . 
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Let zjfa = Z^j (<&) be the double-dimer partition function; for each 
e it is a polynomial function of the entries of the A4. Let 

(6) F($) = lim fV 

~°Z«(Id) 

where Id is the trivial connection. 

Theorem 4. The above limit exists and the function F(<$>) is confor- 
mally invariant, that is, depends only on the conformal type of(U, {z\, . . . 

The proof will show that F($>) is an integral of certain products of 
Green's functions on U. The geometric meaning of F is not clear in 
general; however for the case when U is an annulus see section [9] below. 

Combined with Theorem [3] we have the following corollary. 

Corollary 5. The limiting distribution of the lamination L is confor- 
mally invariant. 

Proof of Theorem [4| We can realize $ by choosing <p m i to be the 
identity on all edges except for the edges crossing a set of disjoint 
simple paths from the Cj to the outer boundary C , as in Figure [2] (in 
blue). For simplicity we choose these edges so that, oriented from the 
boundary to /j, each crossing edge has a white vertex on the left and 
black vertex on the right. We call these collections of edges zippers. 
Let Ei be the set of edges in the ith zipper. 

Let {$t}te[o,i] be a smooth one-parameter family of flat connections 
supported on the zippers with $0 being the identity and $1 = $. It 
suffices to show that 4logZ^($ t ) is a conformally invariant quantity 
plus an error tending to zero as e — > 0. 

Fix e > 0. Let K t be a Kasteleyn matrix associated to <3> t as above, 
with a fixed sign convention independent of t. Let K t be the matrix 
obtained from K t by replacing each SL2(C) entry with its 2 x 2 block 
of complex numbers; thus K t is a matrix of twice the dimension; to 
each vertex v G Q e there are two rows and two columns v^ 1 ' and v ^ of 
K t . Under an infinitesimal change t — > t + 5, Kt+s = K t + 5S t for some 
matrix St supported on the zippers. We have 

|>o g ^»(* t ) = |logQde«-« = ^lo g de t /? t , 
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and 

det(K t + 5S t ) = detK t det(l + 5 StR- 1 ) 

= det K t (l + STiiStR- 1 ) + 0{5 2 ] 



detR 



U + 5J2 Stim^k-^u^m) + 0(6 2 



so that 



^ log det = ^2 S t (u 1 ,u 2 )k t 1 (w 2 ,«i) 



dt 

U\,U2 

Here u±, u 2 run over all rows of k^ . 

Suppose without loss of generality that only A\ changes. Then 
St{ui,u 2 ) = unless U\,u 2 correspond to the vertices of an edge of 
zipper Ei, that is, U\U 2 = w^b^f 1 or U\U 2 = uf'w^ for an edge of 
Ei, and p, q e {1,2}. In this case S t (ui,u 2 ) is equal to K(ui,u 2 ) (the 
edge weight) times either A[(p,q) or (A^ 1 )'^, q) depending on which 
direction u\u 2 crosses E\. 

The sum of the two contributions from an edge = w^bi and its 
reverse is then 

(7) *,(*,„,) Y. d ^K\^H d -^K^W)- 

For simplicity (and without loss of generality, by path invariance) 
let us assume that the path 71 is polygonal with slope ±1, and for 
each segment of 71 the zipper E\ consists of a zig-zag path, alternately 
horizontal and vertical in one of the directions NE,NW,SW,SE (as in 
Figure[2]for the SW path). Suppose moreover that each zig-zag segment 
of Ei has length which is an even number of lattice spacings (except 
perhaps the last segment), that is, an even horizontal and even vertical 
length. 

Let us then compute the contribution to the sum (j7l) for a given seg- 
ment of 71 . Suppose first that it is oriented northeast. Let {(wi, &i)}f=i 
be the corresponding edges of E\. We group these edges into packets 
of four consecutive edges, each packet consisting of two horizontal and 
two vertical edges. 

Consider a horizontal edge wb with white vertex of type Wo, in the 
same packet as a horizontal edge w'b' with white vertex of type W\. 
By Lemma M below, when wb is not within 0(e) of the boundary, 



K-\w, b) = A l {h + eRe(4(^) + F.(z, z)) + 0(e 2 )) 
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and 

K-\w', b') = A x {\l + eHe(4(z) - F_(z, *)) + 0(e 2 )). 

Summing these two contributions gives 

-A 1 + 2A l eReFl(z) + 0(e 2 ). 

For these edges in the reversed orientations the sum is the g-conjugate 
of this (recall that A\ = A^ 1 when det A t = 1): 

l -A- x l +2eReFl*{z)Ai 1 + 0(e 2 ). 
The contribution to the sum ([7]) is then 

{{A^yA^-A^+e ([A^)'A x KeFl{z) + A' 1 ReF|*(z)A r 1 )+0(e 2 ). 
The leading term vanishes: 

leaving the term of order e and a negligible error. 

For the two vertical edges wb and w'b' in the packet we have similarly 

i 



and 



Kf\w, b) = A^I + telm(Fl(z) + F_(z, z)) + 0(e 2 )) 



K-\w', V) = A 1 ( l -I + idm(4(z) - F_(z, *)) + 0(e 2 )). 



Summing these two contributions gives 

% -A 1 + A 1 ielmFl(z) + 0{e 2 ). 

the net contribution to (|7| from these two edges and their reverses is 
then (we must multiply by —i which is the entry K(w,b) = K(w',b'), 
and the leading terms cancel as before) 

e ^(A^yA.lmFliz) + A l \mF ] + \z)A^ + 0(e 2 ). 

Summing the contributions for the packet of four edges we get 



(A^ 1 )' A 1 lm(Fl(z)dz) + A' l hn{F^{z)dz)A^ 1 ) + 0{ 



where we used the notation u dz" to represent e(l + z), the displacement 
from the beginning of the packet to the end of the packet. 
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Now consider the case of a northwest segment of 71. The only change 
is the contribution for the paired vertical edges wb and w'b' '. These are 



K-\w, b) = A 1 {- l -I + ielm(4(*) + F_(z, z)) + 0(e 2 )) 



and 



KrHw', b') = A 1 {~I + <dm(4(z) - F.(z, z)) + 0(e 2 )). 



Summing these and multiplying by i the edge weight and the appro- 
priate matrix we get 

-e ({A^yAJmFliz) + A[lmFl* (z) A^ + 0(e 2 ). 

When added to the horizontal contribution, the net contribution for 
four edges is 

{A?)' A x TMF\(z)dz) + A\\m{Fl*{z)dz)A-[ l ) + 0(e 2 ) 

where u dz" now represents the displacement e(— 1 

In a similar manner the other two possible directions of segments of 
7x also contribute 

(A^ 1 )'A 1 lm(Fl(z)dz) + A' 1 lm(Fl*(z)dz)A^ 1 ) + 0(e 2 ), 

where u dz" represents e times 1 — i, 1 — i according to 

the direction of the segment being NE, NW, SW, SE respectively. 

When wb is within 0(e) of the boundary this formula must be mod- 
ified: the leading terms are no longer of modulus 1/4. However as 
mentioned in the comments after Lemma [8] below, we can choose a lo- 
cal trivialization of the bundle near wb (by isotoping the zipper out of 
the way) to see that the leading terms in K~ 1 (w, b) and w) are 

replaced by the same constant Cf, w and thus still cancel as before. The 
subleading terms of order e only differ from the above when b is within 
o(l) of the boundary, and so when summed these boundary-error terms 
contribute a negligible amount. 

In the limit e — > and the sum becomes the imaginary part of 
contour integral (the 0(e 2 ) term, when summed over the path, is at 
most 0(e) and drops out) 



logdet = Tr j^^Ailm (^j F ] + (z) dz^j + ^Im (^j F ] + *(z) dz^j A{ 1 



d_ 

dt k ^ V</71 / ^ V</71 

These are contour integrals of analytic functions depending only on 
the conformal type of the surface, along the path 71. Thus the limit is 
conformally invariant. □ 
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6. K 1 AND DISCRETE ANALYTICITY 



Let U and U e be as in section 5.1 The goal of this section is to 



determine the asymptotic form of K~ l (b,w) for adjacent vertices b, w 
(Lemma [8] below). In the case of trivial bundle this was worked out in 
[7]. The proof there applies essentially without change to the case of 
a nontrivial flat bundle. We give here an overview of the results of [7J 
and then indicate how the proofs change in the presence of a bundle. 

6.1. Discrete analytic functions. Let V e , E e , F e be the vertices, edges 
and faces of U e , where F e includes the outer face fo and the face fi in- 
side the ith boundary component. A discrete analytic function [2] 
is a function u + iv , where u : V e — > R and v : F e — > IR satisfy the 
discrete Cauchy Riemann equations: for an edge e = X1X2, 

u(x 2 ) - u{x x ) = v(fi) - v(f 2 ), 

where fi, f 2 are the two faces adjacent to edge e, and fi is the face 
to the left when e is traversed from x\ to x 2 - This can be written 
succinctly as 

(9) * du = dv 

where d represents the difference operator and * is the "rotation by 
90°". 

On non-boundary edges this is equivalent to the "discrete Cauchy- 
Riemann equations" 



(10) 



Uy = ~V X 



on, respectively, horizontal and vertical non-boundary edges, where the 
partial derivatives represent discrete differences: u x ((x, y)(x + e,y)) = 
u(x + e, y) — u(x, y) and u y ((x, y)(x, y + e)) = u(x, y + e) — u(x, y). On 
boundary edges these still hold if one interprets v as being constant 
just outside each boundary component. 



A function which satisfies ( 10 ) except at some subset of edges ei, . . . , 
is said to be discrete meromorphic with poles at ei, . . . , e^. In this 
case the defect of the CR equations defines the residue: the residue for 
a horizontal edge x±x 2 where x 2 = x\ + (e, 0) and fi, f 2 are the adjacent 
faces as in ^ is u(x 2 ) —u(x{) —v(fi) +v(f 2 ) and the residue for a ver- 
tical edge xix 2 where x 2 = x\ + (0, e) is i{u(x 2 ) —u(x\) — v(f\) —v(f 2 )). 
Note that the residue is either real (if the edge is horizontal) or pure 
imaginary (if the edge is vertical). Also note that the residue is defined 
even for boundary edges. 
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If u + iv is a discrete analytic function then both u and v are discrete 
harmonic: 

(11) 4«(p) = u(p + (e, 0)) + u(p + (0, e)) + u(p - (e, 0)) + u(p - (0, e)) 

and likewise for v. This follows from summing Q for the four edges 
coming out of a vertex (for u) or the four edges surrounding a face (for 
v). If u + if is meromorphic with a pole at e of residue c G R then u 
is not harmonic at the vertices of e and v is not harmonic at the faces 
adjacent to e. The Laplacian of u at the vertices of e is ±c, depending 
on whether the vertex is the right or left endpoint of e. The Laplacian 
of v is ±c at the upper, resp. lower face. Similar equations hold for 
imaginary residues (at vertical edges). 

If U is a multiply-connected planar domain (or Riemann surface), a 
discrete analytic section of a flat bundle on U e is a section which 
is locally a discrete analytic function in any local trivialization of the 
bundle. 



6.2. Kasteleyn matrix. Recall the definition of the bipartite graph 



Q e from section [5TT] The black vertices of Q e are B = Bq U Bi, where 
Bq are vertices of U e and B\ are faces of U e . The white vertices are 
W = Wo U W\ where Wq are horizontal edges of U e and W\ are vertical 
edges. 

Let K e be the Kasteleyn matrix for Q e whose rows index white ver- 
tices and columns index black vertices, with K e (w,b) = — 1, — i 
according to whether b is adjacent and E, N, W, or S of w. 

If we avoid the boundary, then K € acting on functions on B is the 
discrete d s operator in the sense that u + iv is discrete analytic function 
on U t if and only if u + iv, considered as a function on B (that is, u on 
Bq and iv on B\) is in the kernel of K € . Taking the boundary values 
into account, a function u + iv in the kernel of K e must also satisfy 
v = on the large faces /j and the outer face, and satisfy u = at xq, 
but can have poles at the e^. 

More generally, suppose that (Q e , U) is equipped with a flat bundle, 
and K e the associated Kasteleyn matrix (whose entries are obtained 
by multiplying the above weights by the parallel transports <f> e ). Then, 
away from the boundary, a section is discrete analytic if and only if it 
is in the kernel of K e . 

Lemma 6. As a function of b, if~ 1 (6, w) is a discrete meromorphic 
(matrix-valued) section with poles at the Cj and at w, and zeros at 
b = Xo and b = fj for all j . The pole at w has residue I or il according 
to w G Wq or w G W\ . There is a unique section with these properties. 
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Proof. Fix w; the equations K e (w', b)K~ l {b, w) = 5 WyW iI are linear 
equations for K~ l (b, w), one for each w'. At w' ^ w they are the 
discrete CR equations. For w' a boundary edge or w' adjacent to Xq, 
they correspond to the CR equations if we extend if" 1 (5, w) to be zero 
at b = fi and b = x . The condition on the residue of the pole at w 
is determined by ^ h K e {w,b)K~ 1 {b 1 w) = I. The uniqueness follows 
from invertibility of K t . □ 

6.3. Green's function. The function K~ for the trivial line bundle 
on Q t is related to the Green's function G of the standard Laplacian 
on U e and the Greens function G* on the dual graph U* as follows. 

Lemma 7 ([7], Lemma 9). We have 

G(w + |,6) - G(w - f,6) weW ,beB 
, . , -i(G*(w + *i,b)-G*(w-%i,b)) weW ,beBx 

-i(G(w + |i, 6) - GO - fz, 6)) u> G Wi, b £ B 

where G, G* is the Greens function for U e , U* respectively with respec- 
tively Neuman, Dirichlet boundary conditions. 

With the appropriate definitions of discrete derivatives and 
we can rewrite this as 



dw x dw y ' 



K:\b,w) 



^ weW ,beB 



In [7] the asymptotics of if -1 is written in terms of the continuous 
Green's function, as follows. Let g(u,v) be the analytic function of v 
whose real part is the Dirichlet Green's function g(u,v). Let g*(u,v) 
be the analytic function of v whose real part is the Neumann Green's 
function g*(u, v). Defindj 



and 



dg(u,v) lfdg(u,v) .dg(u,v) 
F + (u,v) = — = - — %- 



du 2 V du x du 



y 



dg(u,v) _ 1 ( dg(u,v) .dg(u,v) 



du 2 V du x du y 



1 Note that this differs from the definition in [7. by a factor of 4; a factor of 2 
is conventional and one is due to the difference in choice of coordinates: here the 
lattice step for the dimer model is e/2, not e. 



18 



RICHARD KENYON 



These are analytic functions of u, v and u, v respectively. It is not hard 
to showBthat 

dg*(u,v) dg*(u,v) 

— (u.v) - 



aw aw 
that is, the difference is only a sign from g. Then for w, b close to w, t> 
respectively with u ^ v we have up to errors of order 0(e 2 ) 

eRe(F + (u,v) +F_(u,v)) weW ,beB 

K-Mh ™\ = I eilm ( F +( u , v ) + u )) w G Wo, b G Bi 

' 1 ' ^ eRe(F+(w,f;) -F_(u,v)) iceWi.fteBi ' 
eilm(F + (u,v) - F-(u,v)) weW 1 ,beB 

Note that for fixed w G Wo, &) is indeed analytic as a function 

of v , except at v = u. Similarly for w G W\. 

In the case w, b are within 0(e) of each other (for example adjacent), 
we define 

fUu) = lim (f+(u,v) -r- — v 

+ v ; v->u\ +K ' ; 2tt(v-u) 

Then 
(12) 

{eRe(Fl(u) + F_ («,«)) u>gWo,&gB 

eiIm(4(w) + F_ («,«)) wGWo^GSi 

6Re(F|( M )-F_Kn)) we^iGfl! ' 

eTm(F|(u) - F_(u,u)) u> G Wi,6 G B 

Here note that is analytic but F_(u,u) is not in general. 

These formulas were shown to hold in the case of trivial bundle but 
the proof applies in the case of a flat SL^C) connection as well; the 
relevant Green's functions are the Green's function of the Laplacian for 
the connection see [6]; it is defined by 

where the sum is over nearest neighbors v' of v. The Green's functions 
G and G* are the inverse Laplacian operators on U e or U* respec- 
tively with the appropriate boundary conditions. For generic $ these 
Laplacians are invertible (see P, Theorem 9] for an expression for their 
determinants) and so G and G* are well-defined. 



2 For the upper half plane g(u,v) = — 5^ log and g*(u,v) — — ^ log(u — 
v)(u — v). For a general simply connection domain U with Ricmann map <f> to U we 

have g(u,v) = -^log |pE^) and 9*i u > v ) = ~h log(0(") _ 0( u ))(0( u ) - 
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Since the Green's function Go on Z 2 satisfies Gq((0,0),v) = — \ for 



any neighbor v of (0, 0) and G ((0, 0), (0, 0)) = [T7j, by g2p we have 
the following result. 

Lemma 8. For an edge wb with w G Wo, crossing the zipper E\ at 
location z G U we have, up to 0(e) 2 additive error, 



A X {±\I + eRe(F ] + (z) + F_(z, z))) weW , b = w±e 

Ail=f{l + ieIm(Fl(z)+F_(z,z))) w G W Q , b = w±i 

Ax{±\l + eRe(Fl(z) - F_(z, z))) w G W x , b = w ± e 

Ai(t{I + ieIm{Fl(z) - F_(z,z))) w G W\, b = w±i 



When z is near the boundary of U (within 0(e)) we have a weaker 
estimate. We first isotope the zipper so that it terminates at some 
other point not close to z\ then the local behavior of K~ x near z is 
that K~ x (b, w) = Ch w + 0(e) for a constant Cb w depending only on 
the distance of edge bw to the boundary and its orientation, see [7j. 
Now if we move the zipper back so that it passes through bw then 
K- 1 (b,w) = A 1 C bw + 0(e). 



7. Peripheral curves 

Theorem [4] applies to multiply connected domains with piecewise 
smooth boundaries. If we are interested in simply-connected domains 
with punctures, one can argue as follows. A peripheral curve is one 
which surrounds a single boundary component or puncture. 

On a multiply connected domain, the above proof shows that for any 
simple curve the probability that it occurs as a curve in the lamination 
of a random double-dimer cover is conformally invariant. This applies 
to both peripheral and nonperipheral curves. In the case when the 
boundary components shrink to points, the number of peripheral curves 
surrounding them tends to infinity. However the nonperipheral curves 
in a lamination are bounded (almost surely) in number and topological 
complexity. Thus one can make sense of the limiting probability of a 
nonperipheral curve. 

The proof of Theorem [1] above shows that the computation of F(<§>) 
only depends on the Green's function of the domain U. In particular 
if we introduce a finite number of holes in U, with small diameter < 5, 
which are not close to any existing boundary or zipper then along any 
zipper the Green's function G changes by an amount tending to zero 
with 6. Since is an integral of G, F($>) also changes by an amount 
tending to zero with 5. Thus the probabilities of homotopy classes of 
double-dimer loops in U are close to those in U' (taking into account 
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only the homotopy classes relative to the boundary of U, not the new 
holes). 

That is, the double dimer loops are insensitive to the addition of 
small holes in the domain (except of course for those small loops which 
come close to or intersect the removed holes). In particular if one wants 
to measure the probability of a homotopy class of nonperipheral loop 
in the complement of a finite number of points, one can simply remove 
small disks of radius 5 around those points and compute the probability 
in the limit as 5 — > 0. 



8. Other graphs 

Theorem [4] applies to other graphs U e as well, on condition that they 
conformally approximate U. 



Let U be as in section |5.1| and let U e be a sequence of graphs which 
conformally approximate U as e — > 0, in the sense that the mesh size 
(diameter of the largest face) tends to zero, and discrete harmonic func- 
tions on U e converge to continuous harmonic functions on U. Equiva- 
lently, the simple random walk on U e converges to the (time-rescaled) 
Brownian motion on U. 

From such a graph U e one can define a bipartite graph G e as above, 
whose white vertices are the edges of U e and black vertices are the 
vertices and faces of U e , see [12]. The Kasteleyn matrix is the adjacency 
matrix with rows indexing the white vertices and columns indexing the 
black vertices, and (in general, complex-valued) signs chosen so that 
around each face (which is a quadrilateral) the alternating product of 
signs is —1, that is, if the face is a quadrilateral ABCD with signs 
a(AB),...,a(DA)),then^0 ) = -l. 

One natural way to choose the signs is to make them complex num- 
bers depending on the geometry of the embedding of Q £ , as follows. 
For an edge wb where b is a vertex of U let K(w,b) = jfz^n that is, 
the unit modulus complex number in the direction of the edge wb. If 
b is a face of U let K(w, b) be the unit complex number perpendicular 
to the edge w and in the direction from the edge to the face b. This 
leads to the Kasteleyn weighting criterion that the four weights of a 
quadrilateral face a, b, c, d satisfy ac/bd = — 1. 

The notion of discrete analytic function on Q e given by (|9| (but not 



(10)) extends to this more general setting, as does the harmonicity of 



the real and imaginary parts (equation (11), with the usual combina- 
torial Laplacian on U). Lemmas [6] and |7| have simple analogs in this 
more general setting: i^ _1 (6, w) is discrete meromorphic with a pole 
of residue e ld at w, when the edge is oriented in the direction 9 (the 
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residue is naturally a 1-form, so depends on a choice of orientation), 
and poles at any other removed edges of U e . In terms of the Green's 
function, when b is a vertex of U we have 

K ib.w) = — aw 

ow 

where again -Jj^ represents a difference operator on U e . A similar for- 
mula with the conjugate Green's function holds when b is a face of 
U. 

In Lemma [8] we have the same formula except that the constants 
1/4 and i/4 are replaced by edge-dependent quantities (but they still 
cancel in ^7|)). The sub leading terms are the same by our hypothesis 
that U t conformally approximates U. We are led to exactly the same 
integral ^ as before. 

9. Annulus 

The case of a cylindrical annulus is particularly easy to compute 
because we can explicitly diagonalize the relevant Kasteleyn matrices. 
Let Q be the graph obtained from a rectangle [0, 2n] x [l,m] C Z 2 by 
identifying for each j vertex (0, j) with vertex (2n, j). This is a bipartite 
graph with 2nm vertices. We use a flat connection with monodromy 
M around the circumference; thus a double-dimer configuration with k 
noncontractible cycles will have weight proportional to (TrM) fe . Since 
the weight only depends on TrM, there is no loss of generality in taking 

M to be a diagonal matrix M — ( ^ \ -i ) ■ Thus we need only consider 



a line bundle (see section 3.3). 

Let a satisfy a 2n = A. We put parallel transport a on horizontal edges 
(x,y)(x + l,y) and 1 on vertical edges. To make a Kasteleyn matrix 
K multiply all vertical edge weights with i = y— 1. If we assume n is 
odd then there are no additional signs needed. 

The eigenvectors of K are then of the following form. Let z, w satisfy 
z 2n = 1 and w 2m+2 = 1. Then for vertices (x, y) G [0, 2n] x [l,m] the 
function 

fz, w ( x ,y) = zX ( wV - w ~ v ) 

is an eigenvector with eigenvalue 

1 ■/ 1 , 
az H h ivw -\ — J. 

az w 

Letting z run over 2n-th roots of unity and w run over the 2m + 2-th 
roots of 1 with positive imaginary part (there are m of these) we have 
all 2nm independent eigenvectors. 
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The determinant of K is then, letting w = e nlk /( m + 1 ) for k 



detK = n ITfa^H h2icos 

11 11 V az m + 1 J 

(az - a k )(az - /3 k ) 



z 2n =l k=l 
m 

nn az 

z 2n =l k=l 



m 

n 

k=l 



(A-0(A-/3f) 



where a k , (3k = i(—cos9 ± y/l + cos 2 9) with 9 = ^tj- Let a k be the 
smaller (in modulus) of the two roots. We can write this as 



detK= (il-PH (n^ + ^X^KI 2 ") 



since a^Pk — 1 and n is odd. Since \a k \ < 1, the terms in the second 
product are negligible for large n unless \a\ ~ 1, that is, except when 
9 « tt/2. If = | + e for small e then \a\ 2n = e - 2n ^ + °^ nt3 \ Suppose m 
is even (the case m odd is similar, see below). Take k = y + j so that 



^ = -JE* = £ + |pt±2 an d UJ2n = e -^^+o(ni/- 3 )_ Xaking 

m+l 2 2(m+l) ' I «■ I o 

large with n/m = r fixed and g = e _T7r we have 



detK = A^JJ(l + gl 2j - 1 lA)(l + g |2j - 1| A- 1 ) 



oo 



n Y[ (i + q 2j + q 3 xf 



3=1 
j odd 

where X = A + A -1 = Tr(M) and where iV is a normalizing factor 
independent of A. 

This gives the probability generating function for the number of loops 
to be (when m is even) 

fpr(»oops)x^ n(;y* +fy. 

S i-i (i + ^ + 'J 2 -') 2 

j odd 

See Figure [3j 
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5 10 15 20 25 30 

Figure 3. Probabilities of to 7 loops as a function of 
1/r = m/n, the ratio of height to circumference of the 
cylinder, in the case m even. 



Similarly one finds, in the case m is odd, the probability generating 
function to be 



£Pr(fcloops)X fc = ^ J] 71 

fc=0 j=2 ^ 

j even 



1 + q j X + q 2 i) 2 



+ q j + q 2j ) 2 



10. Loops surrounding two points 

We show here how to compute the distribution for the number of 
loops surrounding two faces in the upper half plane grid Q = Z 2 fl 
{y > 0}. Let /i,/2 be two faces, and take a flat SL 2 (C)-bundle with 
monodromy A G SL2(C) on a zipper from the boundary to f\ and 
B e SL 2 (C) on a zipper from the boundary to / 2 as in Figure |4j 



-0-3*— 



-g*6 



O-g*- 

— 



FlGURE 4. 
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We let A = jj and B 

2 + e 2 . 

We then have 

oo 2 

Z = QdetK = J2C k (l + € -) k , 

k=0 

where C k is the weighted sum (weighted by 2 C , c being the total number 
of loops) of configurations with k loops surrounding both f\ and f 2 . 
(Note that loops surrounding only one of fi or f 2 are counted correctly 
since TyA = 2 = TrB.) 

The coefficient of e 2 in the expansion around e = of Z is then | of 
the expected number of loops. Higher coefficients give higher moments. 

Let us compute the coefficient of e 2 . We have K = 

where M is the matrix with SL^C) entries and with rows indexing 
white vertices and columns indexing black vertices. Let Mq be the 
corresponding matrix for the trivial connection; the double M is a 
direct sum of two copies of K , the standard scalar Kasteleyn matrix 

for the upper half plane. Now QdetK = V det K = det M and so 
Z = det M = det(M + eS) = det M + e 2 i^o&i' + 

Here u>i&i is an edge of zipper A and w 2 b 2 is an edge of zipper b. (The 
terms of first order in e are zero, as are the terms of order e 2 with both 
edges from the same zipper.) 

However M is just two copies of K , so 

(Mo)5J*g = (detMo)^" 1 ^! 1 ^)^" 1 ^)- 

In the scaling limit for the upper half plane, let w\b\ be a hori- 
zontal edge near a point with complex coordinate z\ and w 2 b 2 be a 
horizontal edge near a point with complex coordinate z 2 . Taking into 
account the orientations (Right or Left) of w\b\ and w 2 b 2 we have 
K,\w\b\)K,\wlbl) = 

if {w\b\,wlbl) is RR 
if {w\b\,wlb\) is RL 
if {w\bl,w 2 2 b\) is LR 
if {w\b\, w 2 bf) is LL 



J J ) for small e. Then Tr(AB) = 



f M\ 
1 M* J ' 



f 4Re(^ + -^)Re(^- + -±r-) 

7T^ ^22—21 22— 21' — 22 21 — 22 ' 

-4lm(^- + -J^)ImM M 

J TT Z ^2 2 — 21 2 2 — Z\> V 21-22 2l~2 2 ' 

-4lm(— MM^" + -^) 

■K 2 - ^ Z2 — 21 22 — 21' ^21 — 22 21—22' 

4Re( — ^)Re(^ M 

V 7T^ ^22—21 22—21 ' ^21—22 21—22 ' 
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Summing over the four possibilities, the sum is a Riemann sum for 
the integral 



detMn 



7T" 



Re 



71 J 12 



1 1 

+ 



(zi - z 2 ) 2 Oi - z 2 ) 2 



dy 1 dy 2 



7l z 



log 



z\ - z 2 



z\ - z 2 



up to errors going to zero with e. 
Thus we have 



Theorem 9. For the double- dimer model on eL 2 P\{y > 0}, the expected 
number of loops surrounding both of the points Z\, z 2 converges as e — > 
to 



4 . 
^log 



7T- 



Z\ - Z 2 



z\ - z 2 



This result can also be obtained using the results of |7J and fact 
that the double-dimer loops are the contours of the height function 
difference of two independent uniform dimer covers. 

An analogous computation for a chordal path can be made: take 
b, w two vertices on the boundary of the upper half plane grid Q = 
eZ 2 n {y > 0} with b < w. Add to Q an edge e connecting b to w. Take 
a double dimer cover of Q, and condition on the event that edge e is 
part of a nontrivial loop. This loop will be a "chordal" double-dimer 
path from b to w in the upper half plane. 

Theorem 10. In the scaling limit on the upper half plane, the proba- 
bility that the chordal double-dimer path from b to w with b < w G R 
passes left of a point z is the harmonic function of z with boundary 
values 1 between b and w and zero elsewhere. 
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